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Centrifugal compressors are compatible with the low exit corrected flows found in the 
high pressure compressor of turboshaft engines and may play an increasing role in turbofan 
engines as engine overall pressure ratios increase. Centrifugal compressor stages are 
difficult to model accurately with RANS CFD solvers. A computational study of the CC3 
centrifugal impeller in its vaneless diffuser configuration was undertaken as part of an effort 
to understand potential causes of RANS CFD mis-prediction in these types of geometries. 
Three steady, periodic cases of the impeller and diffuser were modeled using the TURBO 
Parallel Version 4 code: 1) a k-£ turbulence model computation on a 6.8 million point grid 
using wall functions, 2) a k-£ turbulence model computation on a 14 million point grid 
integrating to the wall, and 3) a k-co turbulence model computation on the 14 million point 
grid integrating to the wall. It was found that all three cases compared favorably to data 
from inlet to impeller trailing edge, but the k-£ and k-co computations had disparate results 
beyond the trailing edge and into the vaneless diffuser. A large region of reversed flow was 
observed in the k-£ computations which extended from 70% to 100% span at the exit rating 
plane, whereas the k-co computation had reversed flow from 95% to 100% span. Compared 
to experimental data at near-peak-efficiency, the reversed flow region in the k-£ case 
resulted in an under-prediction in adiabatic efficiency of 8.3 points, whereas the k-co case 
was 1.2 points lower in efficiency. 
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I. Introduction 

C entrifugal compressors are frequently used in the 
high pressure compressor of turboshaft engines 
due to the low exit corrected flow rates in these 
machines. With the decreasing size of the high 
pressure compressor in turbofan engines as engine 
overall pressure increases, exit corrected flow rates 
may be sufficiently low such that the use of a 
centrifugal stage at the compressor exit may provide 
performance benefits over a typical axial compressor 
stage. 1 Thus, the ability to accurately simulate the 
flow within centrifugal stages is becoming 
increasingly important for compressor designers and 
researchers. However, Reynolds-averaged Navier- 
Stokes (RANS)-based computational fluid dynamics 
(CFD) results published in the literature tend to mis- 
predict the performance of centrifugal compressors 
when compared to available data. 2 * " 5 In general, the 
computations tend to over-predict total pressure ratio, 
efficiency, and in the case of vaned diffuser 
configurations, the choking flow rate. 

The current research effort was undertaken in an attempt to investigate and understand the potential causes of the 
mis-prediction of centrifugal compressor impeller performance by RANS solvers. Figure 1 shows the impeller 
studied in this investigation, a centrifugal compressor known as CC3. This test compressor is a 4.54 kg/s (10 lb m /s) 
machine which was scaled up from a 1.658 kg/s (3.655 lb m /s) Allison Engine design. The design exit corrected flow 
rate of the test compressor is 1.36 kg/s (3.0 lb m /s). The complete design description and hot shape geometry of the 
test compressor are given by McKain and Holbrook. 6 The test compressor was used in aerodynamics studies at 
NASA Glenn Research Center. 7 ' 9 The impeller was tested in both vaned diffuser and vaneless diffuser 
configurations, and detailed measurements of the flow field were obtained. The data available for this test 
compressor qualifies it as a good candidate for further computational studies. 

The CFD solver used in the current work is the TURBO Parallel Version 4 code. TURBO is a three-dimensional, 
viscous, unsteady RANS solver for turbomachinery flows. Flux vector splitting is used in evaluating flux Jacobians 
on the left hand side, while Roe’s flux difference splitting is used to form a higher order TVD scheme for convective 
fluxes on the right hand side. The code can use multiple structured blocks with matching boundaries. Local time 
stepping can be employed to accelerate convergence for steady-state flows. Newton sub-iterations are applied to 

unsteady flows to obtain a converged solution within each 
time step. 10 Multiple blade rows can be modeled with 
sector periodic blade row meshes using a simple periodic 
boundary condition or with a single blade passage per 
blade row using phase lag boundary conditions. 11,12 The 
default turbulence model in TURBO is the NASA/CMOTT 
k-s turbulence model with wall functions. 13 The code has 
previously been validated on various geometries. 14,15 

As a precursor to the current research effort, unsteady 
phase-lag simulations of the CC3 impeller with the vaned 
diffuser configuration were generated using TURBO. The 
resulting speed line at the design speed is shown in Fig. 2 
alongside a comparison to the data. An over-prediction in 
choking flow rate was observed in the CFD result, and this 
speed line serves to illustrate a trend seen throughout the 
literature. 2 " 5 

The CFD results shown are converged, unsteady phase- 
lag computations which were time-averaged over one 
revolution. Computed static pressures at the exit rating 
station were area-averaged. In these simulations, a grid of 
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Figure 2. Computed and experimental total-to- 
static pressure ratios of CC3 vaned diffuser 
configuration at design speed. 




4.3 million points with resolution deemed to be typical of turbomachinery CFD was used. A row of three support 
struts upstream of the impeller, as well as leakage flows, were neglected in the computation. The computation over- 
predicted pressure ratio by 3.3% at the design point inlet corrected flow rate (shown as a filled symbol) and over- 
predicted choking flow rate by 1.9%. The computation was run with ambient inlet pressure of 101,325 Pa, whereas 
the data was taken at approximately 86,000 Pa. The computation was performed with the measured clearances, 
which are minimum tip clearances. If clearances are highly non-axisymmetric in the test article, modeling the 
minimum tip clearance may explain part of the over-prediction in the computed performance. The questions of 
clearance sensitivity and Reynolds number sensitivity are left for future work. Instead, it was hypothesized that the 
over-prediction of performance may be due to inaccurate simulation of the unsteady, non-uniform impeller 
discharge flow and the subsequent impact on associated time-mean aerodynamic blockage at the diffuser throat. 

The scope of the research effort was narrowed to more closely examine the impeller exit flow. The focus of this 
paper is the CC3 impeller in its vaneless diffuser configuration. The compressor was modeled as a steady, periodic, 
single blade row computation in TURBO, which reduced computational expense and run-times. While the vaneless 
configuration obviates the impact of periodic unsteadiness, it was anticipated that lessons learned in modeling this 
simpler case may be useful in future work involving impellers and vaned diffuser configurations. The details of the 
compressor experiments are provided in the next section, followed by details of the computational effort. 


II. Details of the Experiment 



Figure 3. Cross-sectional drawing of the CC3 rig in 
the vaneless diffuser configuration. 7 A is the exit 
rating station, R/R TE =1.18; B is the impeller main 
blade leading edge; C is the impeller splitter blade 
leading edge; D is the leading edge of a row of 
support struts; E is the impeller trailing edge. 
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Figure 4. CC3 inlet survey data measured 2.54 cm 
(1.0 inch) upstream of impeller leading edge. 


Table 1: CC3 impeller design parameters. 


Impeller main blade count 

15 

Impeller splitter blade count 

15 

Design corrected speed 

21,789 RPM 

Design inlet corrected flow rate 

4.54 kg/s 

Design exit corrected flow rate 

1.36 kg/s 

Impeller exit tip speed 

492 m/s 

Inlet blade height 

64 mm 

Exit blade height 

17 mm 

Leading edge tip radius 

105 mm 

Trailing edge radius, R TE 

215.5 mm 

Impeller backsweep 

50° 


The focus of the current work was the CC3 impeller 
in its vaneless diffuser configuration. A cross-section of 
the compressor rig is shown in Fig. 3. Table 1 
summarizes important design parameters. The impeller 
consists of 1 5 main blades and 1 5 splitter blades having 
50 degrees of backsweep from the radial direction. The 
design corrected speed was 21,789 RPM, and the 
measured minimum clearances at this speed were 
0.1524 mm (0.006 inches) at the leading edge, 0.6096 
mm (0.024 inches) at the midchord “knee” and 0.2032 
mm (0.008 inches) at the trailing edge. 7 At this 
clearance schedule and in the cold shape, the impeller 
hub at the trailing edge sat axially forward of the 
diffuser hub by 0.127±0.0508 mm (0.005±0.002 
inches). This step in the hub was neglected in the 
computations. 

Total pressure and total temperature exit rakes were 
located at eight evenly spaced circumferential locations 
at radius ratio R/R TE =1.18. Each total pressure rake had 
four elements located at 19.2%, 43.4%, 68.2%, and 
93.1% of span, and each total temperature rake had 
three elements located at 23.1%, 55.9%, and 88.9% of 
span. Performance parameters were calculated using 
arithmetic averages of these rake elements. An inlet 
probe was used to acquire spanwise total pressure at a 
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station 2.54 cm (1.0 inch) upstream of the impeller main blade leading edge (Fig. 4). Spanwise total pressure 
surveys were also acquired by using a constant-blockage probe immersed downstream of the impeller trailing edge 
at radius ratio R/R TE =1.10. Features which were neglected in the computations include a row of three inlet struts 
upstream of the impeller, impeller blade fillets, and backplate leakage paths. A set of 16 equally spaced 5/1 6 th inch 
diameter pins, which held together the diffuser hub and shroud, were located in the test compressor at radius ratio 
R/R te = 1.45. These structural elements were thought to be sufficiently far downstream of the exit rating station and 
were neglected in the computational model. 


III. Details of the Computations 

The default turbulence model in TURBO Parallel 
Version 4 is the NASA/CMOTT k-s model with wall 
functions. 13 In addition to the default k-s model, 
Wilcox’s k-co turbulence model was used. 16 The 
implementation of the k-co model does not include wall 
functions, and the k-co model required a grid with 
sufficiently low y + so as to allow integration to the wall. 
Thus, two different grids were used in the current effort. 
A “baseline” grid with approximately 6.8 million points 
and y + on the order of 50 was used for k-s computations 
with wall functions enabled. A “refined” grid containing 
approximately 14 million points was also generated, 
allowing integration to the wall for k-co as well as k-s 
_ „ TT . + , computations, in an effort to characterize the effect of 

igure . istogram comparing y o -s an -co var y^ n g g r ^ resolution on aerodynamic performance. A 
computations on the refined grid near design point. histogram of y - for the k _ 8 and k _ 0) computations on the 

refined grid is shown in Fig. 5. Typically, y + was less than 3 for a majority of the walls. Values of y + approached 5 
along the blunt trailing edge of the impeller and at the shroud near the trailing edge of the suction side of the 
impeller, but these were within the viscous sublayer region. In all computational cases, inlet turbulence intensity was 
set to 1%, and inlet turbulent eddy viscosity ratio was set to 100.0. These values correspond to an inlet turbulence 
length scale of approximately 0.6% of the impeller inlet blade height. 

TURBO solves the unsteady RANS equations on multiblock structured grids. All grids, boundary condition files, 
and block connectivity files for TURBO were generated using the Turbomachinery Gridding System (TGS). 17 The 
measured tip clearances, which were minimum tip clearances, were included in the computational grid. Impeller 
blades and tip clearances were gridded with O-blocks, and H-blocks were used to connect the O-blocks around the 
impeller main and splitter blades. The block topology of the refined grid is shown in Fig. 6. The computational 
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Figure 6. Block topology of the CC3 impeller and 
vaneless diffuser. 



Figure 7. Meridional view of the CC3 impeller and 
vaneless diffuser model. 
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domain consists of the axial inlet, an impeller and splitter passage, the vaneless diffuser, and the S-bend to the axial 
outlet. A meridional view of the computed flow path is shown in Fig. 7. 

Steady, single passage periodic simulations were converged at design speed using local time stepping. Solid 
walls were modeled with no-slip boundary conditions and were assumed to be adiabatic. The specific heat 
coefficient ratio was assumed constant. The inlet boundary condition was specified with spanwise profiles of total 
pressure, total temperature, radial flow angle, and tangential flow angle. The total pressure profile, shown in Fig. 4, 
was the measured inlet probe survey data from the CC3 experiment. The computations were referenced to an inlet 
pressure of 101,325 Pa, whereas the data was taken at reduced inlet pressures on the order of 85,500 Pa. Inlet total 
temperature was assumed constant 288.15 Kelvin, and zero inlet tangential flow angle was assumed. The profile of 
inlet radial flow angle was found by linear interpolation between the geometry angles of the hub and shroud at the 
inlet. The computations were throttled from choke, up the speed line, but were not taken to the stall boundary. The 
operating point was set by using a radial equilibrium exit boundary condition, specifying exit static pressure at the 
shroud. 


IV. Results 

Design speed computations of the impeller and vaneless diffuser were generated using TURBO. Three computed 
speed lines are compared to data in Fig. 8. The three computed results are the k-s result with wall functions on the 
baseline grid, and the k-s and k-co results integrating to the wall on the refined grid. The computations were not 
taken to the stall boundary or deep into choke, as comparisons to the near-peak-efficiency data point were of interest 
in the current effort. The computed results were area-averaged in the circumferential direction at the exit rating 
station location, R/R TE =1.18. The area- averaging of the steady computations in the rotating frame of reference is 
analogous to the time-averaging done by the steady rake instrumentation, which is fixed in the stationary frame of 
reference while immersed in the unsteady rotor exit flow. The circumferential averaging yielded a spanwise profile 
of computed total pressure and total temperature. An arithmetic average of the values of total pressure and total 
temperature at those spanwise locations which correspond to the rake element locations was calculated, which was 
analogous to the experimental averaging approach. 

The choking flow rates of the baseline and refined k-s computations agree to within 0.33% of each other. The k- 
cd computation choked at corrected flow 0.5% lower than the refined k-s computation. This difference is attributed 
to a slightly larger growth in the boundary layer in the k-co result along the pressure side of the main impeller blade. 
The larger boundary layer results in increased aerodynamic blockage, reducing both effective flow area and choking 
flow rate. The k-co computation under-predicts choking flow rate by 1 .22% as compared to data. 

Both k-s results under-predict total pressure ratio and efficiency as compared to the data. The k-co computations 
better match the data on these performance parameters. The near-peak-efficiency data point and the computed 
operating points at nearby inlet corrected flow rates are shown as filled symbols in Fig. 8. These operating points 
were compared. All three sets of computations over-predicted the work done by the impeller, as shown by the total 
temperature curves which are approximately 0.6% higher than the data. Compared to the near-peak-efficiency data 
point, the baseline k-s result was 7.2 points lower in efficiency and 7% lower in total pressure ratio. The refined k-s 
result was 8.3 points lower in efficiency and 8% lower in total pressure ratio than the data. The k-co result was 1.2 
points lower in efficiency and 0.6% higher in total pressure ratio than the data. In this case, the choice of turbulence 
model evidently had a much greater impact on the flow field than the selected level of grid density. 

A comparison of shroud static pressures at these operating points is shown in Fig. 9. The data have been 
normalized by the plenum pressure of 85,360 Pa and the computations are referenced to inlet total pressure 101,325 
Pa. The shroud static pressures of the data, the k-co case, and both k-s cases show very good agreement from the 
inlet to the impeller trailing edge. Downstream of the trailing edge, the baseline and refined k-s computations have 
virtually identical shroud static pressures. The k-s computations have reduced pressure recovery through the diffuser 
as compared to the data and to the k-co computation. In the vaneless diffuser, the k-co simulation matches the data 
well through the initial half of the diffuser, but over-predicts pressure recovery towards the S-bend and beyond. 

In Fig. 10 spanwise profiles of the computed flow angles (from tangential) and total pressures are compared 
against the profiles measured by the constant-blockage probe at radius ratio R/R TE =1.10. This is a location nearly 
midway between the impeller trailing edge and the exit rating station, where the shroud static pressures of the k-co 
case and the k-s cases are diverging. The uncertainty of the measured flow angle is not obvious, so measured flow 
angles are shown with both ±1° and ±3° error bars for reference. The comparison shows that the baseline and refined 
k-s computations both yield nearly identical total pressure and tangential flow angle profiles. The flow angle profiles 
indicate that the k-s computations have a region of reversed flow (i.e. negative tangential flow angle) in the top 25% 
of the passage. This large region of reversed flow is an aerodynamic blockage which causes higher velocities and 
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Figure 8. Total pressure ratio, total temperature ratio, and adiabatic efficiency of CC3 vaneless diffuser 
configuration at design speed. Note: the computed speed lines were not run to stall. 


thus higher flow angles in the bottom 75% of the passage, as compared to the data. The reversed flow near the 
shroud is evidently the cause of the large deficit in total pressure observed in the k-s computations, which disagree 
with the measured pressure profile above 55% span. The data indicates reversed flow within the top 10-15% of span 
near the shroud. The flow angle profile from the k-co computation agrees with the data in this regard. The k-co 
computation over-predicts flow angle in the bottom 50% of span, but to less of a degree than either k-s computation. 
The k-co computation misses the pressure deficit seen in the survey data above 75% span. 

In Fig. 11, a comparison of flow angles is made between the three computations at two locations; just 
downstream of the trailing edge (R/R TE =1.01), where the shroud static pressures are in agreement, and at the exit 
rating station (R/R TE =1.18), where there is a large difference in shroud static pressures. The flow angle profiles of 
the three computations match well near the trailing edge. At R/R TE =1.01, the flow angle profiles of the k-s cases 
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Figure 9. Computed and experimental static pressure ratio along the shroud surface. Vertical lines mark 
geometry features of interest. 


match each other to within ±0.5°, and the k-co case matches the k-s cases within ±1.5°, with slight disagreement 
above 95% span. The k-s cases indicate reversed flow beginning at 97% span, the k-co case indicates reversed flow 
at 95% span. There is large disparity between the k-s cases and the k-co case at radius ratio R/R TE =1.18. The flow 
angle profiles of the two k-s computations, which agree to within ±1.0°, indicate that the flow has reversed in nearly 
30% of the span near the shroud. The k-co case has reversed flow in only the upper 5% of span near the shroud. 

The deficit in static pressure recovery seen in the k-s computations is attributed to the large region of reversed 
flow. This region is shown in Fig. 12, which plots the computed radial velocity contours near the impeller trailing 
edge, R/R te = 1.01, and at the exit rating station, R/R TE =1.18. The contours near the trailing edge show the splitter 
blade in the middle of the plot, and the main blade at the left and right sides of the plot. At R/R TE =1.01, the low- 
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Figure 10. Area-averaged computed total pressure ratio (left) and flow angle (right) compared to survey data 
at radius ratio R/R T f=1.10. 
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Figure 11. Area-averaged flow angles comparing results of the k-£ and k-co turbulence models at radius ratios 
R/R te = 1.01 (left) and R/R TE =1.18 (right). 

momentum regions near the hub that extend into mid-span are associated with the impeller blade wakes. The low 
radial velocity region at the shroud is the tip clearance flow. The spanwise extent of the clearance flow is slightly 
larger in the k-s computations as compared to the k-co computation, indicating a stronger tip clearance vortex in the 
k-s computations. The flow fields of the two different turbulence models are qualitatively and quantitatively similar 
at the trailing edge. If performance is computed from the inlet plane to the plane defined by R/R TE =1.01, the refined 



Figure 12. Radial velocity contours (normalized by 287.6 m/s) at the trailing edge of the impeller, R/R TE =1.01 
(left), and at the exit rating station R/R TE =1.18 (right). Impeller blade rotation is right-to-left. 
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k-s and k-co computations agree to within 0.7% in total pressure ratio and to within 1.0 point in adiabatic efficiency. 
At the exit rating station R/R TE =1.18, the two k-e computations are qualitatively very different from the k-co 
computation, and the computed performance parameters at R/R TE =1.18 amount to a 8.4% difference in total pressure 
ratio and 7.0 point difference in adiabatic efficiency between the refined k-s and k-co computations. 

Figure 13 plots entropy on planes of constant radius through the last 30% of impeller chord to show loss 
development through the passage. In the figure, the pressure side of the main blade is on the right, the splitter blade 
is in the middle, and the suction side of the main blade is on the left. The computations differ at the corner formed 
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Figure 13. Entropy contours at constant-radius planes over the last 30% of impeller chord, with the splitter 
blade shown in the center of the passage. 


by the pressure side of the main blade and the 
shroud. The higher magnitude of entropy in the 
refined k-s case indicates larger loss due to a 
stronger tip clearance flow than in the baseline 
k-s case. The spanwise extent of the high 
entropy region is smaller in the k-co case, 
indicating a weaker tip clearance flow than in 
the k-s cases. Qualitatively, the entropy 
contours among the three computations are 
similar through the trailing edge. 

Figure 14 plots radial velocity through the 
radial diffuser and axial velocity after the S- 
bend. These velocities were area-averaged in 
the circumferential direction for the k-s and k-co 
results on the refined grid. This shows the 
extent of the reversed flow through the diffuser 
along the shroud in each case. In both 
computations, the impeller exit flow is similar, 
but differences arise in the diffuser. Both 
computations show reversed flow at the shroud 



Figure 14. Area-averaged velocity contours (normalized by 
287.6 m/s) showing streamwise components (radial 
component in the radial diffuser and axial component after 
the radial-to-axial S-bend). 
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at close proximity to the impeller trailing edge. In the k-e computation, the reversed flow persists up the radius of the 
diffuser along the shroud surface. In the k-co case, there is a relatively small bubble of reversed flow which weakens 
as radius increases through the diffuser. The separated region in the k- 8 cases reattaches through the S-bend and 
there is no negative axial velocity at the end of the computational domain. The diffuser flow field had large 
sensitivity upon the choice of turbulence model. A discussion on the differences in the diffuser flow field is 
presented next, as well as comments and suggestions for future work. 


V. Discussion 


A better understanding of the cause of the disparate diffuser flow fields between the k-e and k-co computations is 
of interest for future work. The flow fields indicate that the k-co computation has greater turbulent mixing in the 
diffuser, which causes the shroud separation to weaken and mix out in a manner that is not occurring in the k-e 
computations. Figure 15 shows a band of high turbulence kinetic energy (TKE) in the diffuser of the k-e case which 
is not evident in the k-co case. This region is attributed to the shear layer at approximately 70% span due to the 
interface of the separated flow and core flow in the k-e case. Generally, it was observed in the k-co computation that 
the levels of both production and dissipation of TKE were higher than the k-e computations near the shroud within 
the impeller passage, but these terms were lower than the k-e computations through the diffuser. 



Figure 15. Turbulence kinetic energy at constant-radius planes from R/R TE =0.9 to R/R TE =1.18. The splitter 
blade is shown in the center of the passage. 



Figure 16. of the refined grid k-e computation 
at radius ratio R/R TE =1.18 (exit rating station). 


It may be of interest in future work to generate a 
simulation using a standard k-e model, which sets constant 
C M =0.09 in the turbulent eddy viscosity equation rather than 
a computed value as in the CMOTT k-e model. As shown in 
Fig. 16, the computed in the refined CMOTT k-e result 
yields values lower than 0.09 in the shear layer region at 
70% span in the diffuser at the exit rating station. A standard 
k-e model would increase in this region, resulting in 
increased turbulent eddy viscosity and turbulent mixing, 
which may mix out or weaken the separation observed in the 
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current k-s cases. 

The current results suggest that the flow developing through the radial diffuser may largely depend on the 
selection of turbulence model. In a vaned diffuser configuration, a strong separation entering the diffuser impacts the 
aerodynamic blockage, which affects the pressure rise and flow capacity in the diffuser. The impact of different 
turbulence models on this separation, and thus on the stage performance, should be analyzed in future work. The 
vaned diffuser computation poses a more complicated problem, as the blockage will also be dependent upon the 
unsteadiness due to the impeller/vane interaction in a vaned diffuser configuration. 

VI. Conclusion 

A computational study of the CC3 impeller in its vaneless diffuser configuration was undertaken. TURBO CFD 
simulations were generated on: 1) a baseline grid using the default CMOTT k-s turbulence model, 2) a refined grid 
using the default CMOTT k-s turbulence model, and 3) a refined grid using Wilcox’s k-co turbulence model. The 
results indicated that all three computational cases matched the data well from the inlet, through the impeller, and up 
to the trailing edge. The area-averaged flow angle profiles of the baseline and refined k-s cases matched to within 
±0.5° just downstream of the impeller trailing edge (R/R TE =1.01). At this radius ratio, the flow angle profile of the 
k-co case matched both k-s cases to within ±1.5°, and the refined k-s and k-co computations agreed to within 0.7% in 
total pressure ratio and to within 1.0 point in adiabatic efficiency. Beyond the trailing edge and into the diffuser, the 
k-co computation continued to match the data, but the k-s computations yielded pressure and efficiency deficits as 
compared to data due to a large reversed flow region extending from the shroud surface of the vaneless diffuser. At 
the experimental exit rating station (R/R TE =1.18), the baseline k-s computation is 7.2 points lower in efficiency and 
7% lower in total pressure ratio than the data. The refined k-s computation is 8.3 points lower in efficiency and 8% 
lower in total pressure ratio than the data. The k-co computation is 1.2 points lower in efficiency and 0.6% higher in 
total pressure ratio than the data. It was shown that the potential impact of the choice of turbulence model may be 
quite large on the diffuser flow field. 
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